Striped phase in a quantum XY-model with ring exchange 
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We present quantum Monte Carlo results for a square- lattice 5* = 1/2 XY-model with a standard 
nearest-neighbor coupling J and a four-spin ring exchange term K. Increasing K/ J, we find that the 
ground state spin-stiffness vanishes at a critical point at which a spin gap opens and a striped bond- 
plaquette order emerges. At still higher K/J, this phase becomes unstable and the system develops 
a staggered magnetization. We discuss the quantum phase transitions between these phases. 
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Ring exchange interactions have for a long time been 
known to be present in a variety of quantum many-body 
systems p| and have been investigated rather thoroughly 
in solid 3 Hc ||. They are also important for electrons in 
the Wigner crystal phase || 0|. In strongly correlated 
electron systems, such as the high-T c cuprates and re- 
lated antiferromagnets, ring exchange processes are typ- 
ically much weaker than the pair exchange J and 
are often neglected. Four-spin ring exchange has, how- 
ever, been argued to be responsible for distinct features 
in the magnetic Raman || and optical absorption spec- 
tra |7) . Neutron measurements of the magnon dispersion 
have also become sufficiently accurate to detect devia- 
tions from the standard pair exchange Hamiltonian (the 
Heisenberg model) and such discrepancies have been at- 
tributed to ring exchange || ||] . Recently, ring exchange 
has attracted interest as a potentially important interac- 
tion that could lead to novel quantum states of matter, 
in particular 2D electronic spin liquids with fractional- 
ized excitations 0, |ll|, ||, [d| |||. Furthermore, for 
bosons on a square lattice ring exchange has been shown 
to give rise to a "exciton Bose liquid" phase Q . 

Here we study the effects of ring exchange in one of the 
most basic quantum many-body Hamiltonians — the spin- 
1/2 XY-model on a 2D square lattice. We use a quan- 
tum Monte Carlo method (stochastic series expansion, 
hereafter SSE [p7Hl8[[l9|) to study the low-temperature 
behavior of this system including a four-spin ring term. 
Defining bond and plaquette exchange operators 



Ba = S+S- 



SrS+ = 2(S?S* + SVS?) 



ijkl — &j °l "1" &i Sj^S k Sf, 



Pijki — Si^S, 

the Hamiltonian is 
H = - 



(ij) (tjki) 
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where (ij) denotes a pair of nearest-neighbor sites and 
(ijkl) are sites on the corners of a plaquette. For K = 
this is the standard quantum XY-model, or, equivalently, 



hard-core bosons at half-filling with no interactions apart 
from the single-occupancy constraint. This system un- 
dergoes a Kosterlitz-Thouless transition at T/J m 0.68 
pfj| , pl| and has a T = ferromagnetic moment M x — 

(Sf) 
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23( | . The if -term corresponds to re- 
taining only the purely x- and y-terms of the full cyclic 
exchange. 

In a soft-core version of the pure ring model (J = 0), 
Paramekanti et al. recently found a compressible but 
non-superfluid phase (exciton Bose liquid) for weak on- 
site repulsion U p6| . As the hard-core limit is ap- 
proached they found a transition to a staggered charge- 
density-wave phase. Hence, the ground state of the spin 
Hamiltonian (^|) can be expected to change from an easy- 
plane ferromagnet with a finite spin stiffness p s and a 
magnetization (M x ) at low Kj J to an Ising-like antiferro- 
magnet with vanishing p s and a staggered magnetization 
(M s ) = (- 1)^+2" (Sf) at large K/J. The central result 
of our simulations is that the competing J and K interac- 
tions give rise to yet a third phase at Kj J ~ 10; a striped 
bond-plaquette phase where the expectation values (fly ) 
and (Pijki) alternate in strength with a period of 2 lat- 
tice spacings in one of the lattice directions. An example 
of this order is illustrated in Fig. |l|. A similar columnar 
"bond charge" phase was recently predicted based on a 
lattice field-theory including a plaquette term The 
field theory also has fractionalized phases, of which we 
have found no evidence. Hence, the microscopic mecha- 
nisms leading to fractionalized spin liquids remain to be 
clarified. 

The SSE simulation method |l7], [ll| 19 that we use 
here has previously been applied to a variety of spin and 
boson models with two-particle interactions, including 
the Hamiltonian (§) with K = (the XY-model) ||. 
The generalization to include the four-spin if -term is 
relatively straight-forward, although non-trivial new pro- 
cedures had to be developed for large- Kj J simulations 
p4j . Bond and plaquette strengths such as those shown 
in Fig. [l] were obtained using open-boundary rectangu- 
lar L x x L y lattices with L y = 2L X . The translational 
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FIG. 1: Plaquette (left) and bond (right) strengths at the 
center of a 64 x 128 open-boundary lattice at K/J — 10 and 
T = J/8. The plaquette strengths are represented by shades 
of gray with the weakest (Pyw) = 0.222 (white squares) and 
strongest 0.468 (black squares). The bond strengths are in- 
dicated by the width of the line segments, with the weakest 
(Bij) = 0.181 and the strongest 0.505. 



and rotational symmetries are then broken and a unique 
static bond-plaquette strength pattern can be observed 
when K/J ~ 10 at T/J < 0.5. For K/J < 8 no order 
is visible at the centers of large lattices at any tempera- 
ture. The modulations seen within the stripes in Fig. |l| 
are strongest at the four corners of the lattice and de- 
crease as the center is approached. They also decrease as 
the lattice size is increased and in the thermodynamic the 
striped state should therefore be analogous to the four- 
fold degenerate columnar spin-Peierls state of Ref. [TJ|. 

Our conclusion that the stripes are stable is based 
on finite-size scaling of correlation functions on periodic 
L x L lattices. The striped phase can be detected using 
the bond or plaquette correlations. Here we consider the 
plaquette structure factor 
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where a\, . . . , are the sites belonging to plaquette a. 
We have studied the full q-dependence and only found 
peaks at (0, 7r) and (71", 0). Hence, the modulations within 
the stripes seen in Fig. |l| are indeed induced by open 
boundaries. The spin structure factor is defined as 



(5) 



where = (Xi,yi) is the lattice coordinate. We will 
analyze the staggered and striped order parameters per 
site, defined as 



(M 2 P ) 



S(ir,n)/L 2 , 
P(ir,0)/L 2 . 
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FIG. 2: Spin stiffness and plaquette-stripe order parameter 
vs ring-exchange coupling for a 64 x 64 system with periodic 
boundary conditions at T = J/8. 



The spin stiffness (the superfluid density in the boson 
representation) is defined by 



d 2 E{cp) 



(8) 



where E(<fi) = (H(4>))/L 2 and the twist <j> is imposed 
in the x or y direction so that the corresponding bond 
operators (@) become B^) = cos {4>){Sf S] + SfS]) + 
sin (<j))(SfS$ - SVSf). The derivative at <j> = in Eq. (§) 
can be directly estimated using the winding number fluc- 
tuations in the SSE simulations p5|. 

Fig. |2| shows the spin stiffness and the stripe order pa- 
rameter on an L = 64 lattice at T = J/8 (where the 
results are almost converged to their ground state val- 
ues). The stiffness becomes very small at K/J w 8, 
where the stripe order increases significantly. Finite-size 
scaling shows that the stripe order survives in the ther- 
modynamic limit. Results for K/J = 8.5 and tempera- 
tures sufficiently low to give the ground state are shown 
in Fig. ||. For L > 32 the data graphed versus 1/L fall 
on a straight line, which extrapolates to a non-zero value 
as L — > 00. Based on results |l6| for the soft-core version 
of the J = model (or K — > 00) the staggered mag- 
netization can be expected to be non-zero for large K. 
However, as also shown in Fig. [| at K/J — 8.5 (M|) 
decreases as 1/L 2 for large lattices, implying that the 
spin-spin correlations are short ranged [S(n, tt) is finite]. 
Fig. H also shows results for K/J = 64, where the scaling 
behaviors of the two quantities is reversed — (Mp) decays 
as 1 /L 2 whereas (Mf ) extrapolates to a non-zero value. 
Note that the size-dependence of Mj is non-monotonic, 
with a minimum around L « 10. Such non-monotonicity 
has previously been observed for a spatially anisotropic 
spin model [M where it was attributed to the presence of 
two different low-energy scales in the system. The non- 
monotonicity seen at K/J = 64 in Fig. |^ indicates that 
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FIG. 3: Finite-size scaling of the ground state staggered mag- 
netization (open circles) and the plaquette-stripe order pa- 
rameter (solid circles) at K/J — 8.5 and 64. The dotted 
straight lines show extrapolations of the infinite-size order pa- 
rameters. The dashed curves show the form ~ 1/L 2 expected 
asymptotically when there is no long-range order. 

the stripe correlations remain strong with a correlation 
length ~ 10 lattice spacings. The location of the mini- 
mum in (Mj) moves to lower 1/L as K/J is decreased, 
indicating growing stripe correlations. The strong stripe 
correlations in the staggered phase makes it difficult to 
determine the (Ms) versus K/ J curve. Our simulations 
show that the stripe order persists at least for K/J up 
to 12, and also that the staggered correlations are short 
ranged up to this coupling. The stripe correlations are 
short ranged for K/J > 16. Between K/J — 12 and 16 
the two phases could either co-exist or be separated by a 
first-order transition. Simulations of larger lattices will 
be required in order to clarify the interesting transition 
region. 

The superfluid-striped transition appears to be of sec- 
ond order, although we cannot exclude a very weakly 
first-order transition (which was argued to be more likely 
in Ref. [ll|). The vanishing of the spin stiffness seen in 
Fig. ^ indicates the opening of a spin gap. A spin gap 
can be inferred also from the temperature dependence of 
the uniform magnetic susceptibility, 

*» = i4((l>") 2 )- < 9 » 
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FIG. 4: Temperature dependence of the uniform magnetic 
susceptibility for L = 80 systems close to the superfluid- 
striped transition. Statistical errors are of the order of the 
size of the symbols. The dashed line shows the linear behav- 
ior expected for a quantum phase transition with = 1. 

Fig. ^ shows the T-dependence for L = 80 (sufficiently 
large to eliminate finite-size effects). The T — > suscep- 
tibility vanishes for K/J between 7.90 and 7.95, i.e., a 
spin gap is present above a critical coupling in this range. 
The temperature independence of Xu at K/J = 7.80 at 
the two lowest temperatures is expected on account of 
this being the behavior in the XY- model jl8], ^3l p7| . 
The behavior for K/J = 7.90 and 7.95 is consistent with 
Xu ~ T at the critical coupling, which is indicative of 
a T = quantum critical point with dynamic exponent 

z = 1 |H, EH- At K / J ~ 7 - 9 we have verified that the 
stripe structure factor indeed exhibits non-trivial finite- 
size scaling, P(7T,0) ~ L e with e < 2, but the statistical 
accuracy is not sufficient for determining the exponent to 
a meaningful precision. Nevertheless, power-law scaling 
for the same K/J at which the spin gap opens supports a 
continuous quantum phase transition with no intervening 
disordered phase or co-existence region. 

In summary, the spin- 1/2 XY- model with ring ex- 
change exhibits three different ground state orderings 
as a function of the strength of the ring term. The 
superfluid-striped transition appears to be a continuous 
quantum phase transition, whereas the striped-staggered 
transition most likely is of first order. Since the sign of 
the J-term in (^) is irrelevant (the sign of the if-term 
is relevant) the superfluid-striped transition could pos- 
sibly, in an extended parameter space, connect to the 
order-disorder transition in the two-dimensional Heisen- 
berg antiferromagnet with frustrating interactions |sT| . 
We also note that the staggered-striped-superfiuid phase 
behavior versus J/K shows interesting similarities to the 
high-T c cuprates, where the pseudogap phase interven- 
ing between the antiferromagnetic and superconducting 
phases exhibits strong stripe correlations p(J. Although 
the microscopic physics and symmetries are clearly differ- 
ent, a detailed study of the staggered-striped transition 
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may still be useful in this context. 

In spite of the absence of a spin liquid phase, the pres- 
ence of three distinct ordered ground states, and the 
phase transitions between them, puts the J—K model (||) 
in an important class among the basic quantum many- 
body Hamiltonians. Although the interesting large--ftf/J 
region may not be of direct relevance to real systems, we 
expect this and related model to be very useful as sys- 
tems where complex quantum states and quantum phase 
transitions can be further explored on large lattices with- 
out approximations. Although other models, such as the 
frustrated Ji — J 2 Heisenberg model [31 , may show simi- 



lar or potentially even more complex behavior, sign prob- 
lems affecting quantum Monte Carlo makes it difficult to 
obtain conclusive results. It would clearly be interested 
to study also the J — K model with a positive sign for 
the K-term, in particular to determine whether fraction- 
alized spin liquid phases could arise, but unfortunately 
this also leads to sign problems. 

The J — K model with the sign of K chosen here can 
be modified in several interesting ways and still be eas- 
ily accessible to simulations using the SSE method. For 
example, when relaxing the hard-core constraint there 
should be a transition to a exciton Bose liquid phase jl6) , 
both as a function of on-site repulsion U for large K/J 
and as a function of K/J. It will also be interesting 
to include a magnetic field to "dope" the striped and 
staggered phases. Transitions between different charge- 
density phases and the question of the existence of doped 
supersolid phases have recently been studied numerically 
for boson models where charge-density phases are sta- 
bilized due to diagonal density-density interactions |$2| . 
In contrast, the striped phase found here arises out of a 
competition between two kinetic terms and it may hence 
behave differently upon doping. 
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